2540 km: Bimagic baseline for neutrino oscillation parameters 
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We show that a source-to-detector distance of 2540 km, motivated recently [l] for a narrow band 
superbeam, offers multiple advantages for a low energy neutrino factory with a detector that can 
identify muon charge. At this baseline, for any neutrino hierarchy, the wrong-sign muon signal is 
almost independent of CP violation and 613 in certain energy ranges. This allows the identification 
of the hierarchy in a clean way. In addition, part of the muon spectrum is also sensitive to the CP 
violating phase and 813, so that the same setup can be used to probe these parameters as well. 
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Introduction. — The data from ongoing neutrino ex- 
periments confirm that neutrinos have distinct masses 
TOi,m2,m3 and the three neutrino flavors v^.^v^^Vt mix 
among themselves. While the mass squared difference 
Am^i = rn^—rnl and the magnitude of Ato§]^ = ra^—rnl, 
as well as two of the mixing angles, 9i2 and 6*23, are well 
measured, three parameters of the leptonic mixing ma- 
trix still remain elusive: the mixing angle ^13, the sign 
of Am|]^, and the CP phase Sqv The determina- 
tion of hierarchy (NH/normal: Ato§]^ > 0, IH/inverted: 
Amgj^ < 0), in particular, would be crucial in identifying 
the mechanism of neutrino mass generation [3 . 

If the actual value of ^13 is not much below the current 
3(7 bound of ^13 < 12°, it may be measured at detectors 
at a distance of ^ 1 km from a reactor/ accelerator. In 
order to determine the hierarchy, however, the most effi- 
cient avenue is to have the neutrinos travel through Earth 
for thousands of km before detection. Here, the differ- 
ence between Earth matter effects in the two hierarchies 
can help in distinguishing them. This can be achieved, 
for instance, by using the decay of accelerated muons - 
/i"*" or - as a source ("neutrino factory" (NF)) and 
a detector that can detect muons and identify them as 
right-sign (the same sign as the source) or wrong-sign. 

The wrong-sign muon signal is hailed as the "golden 
channel" since it is sensitive to all the three parameters: 
6*13, the sign of Amgj^, and (5cp- However, the dependence 
on (5cp also introduces large uncertainties, making the 
unambiguous determination of the true parameters diffi- 
cult [U Ej . A potential way out is to have the detector at 
~ 7500 km ("the magic baseline" 01^) from the source, 
where the effect of CP violation vanishes for both hierar- 
chies. However, this very feature makes it impossible to 
measure the CP phase at this baseline. Moreover, such a 
long baseline requires an extremely well-collimated muon 
source, else the flux at the detector is highly reduced. 

It is therefore desirable to look for a shorter baseline 
that will still give a wrong-sign muon signal independent 
of the CP phase for one of the hierarchies, albeit only 
in a part of the spectrum. The remaining part of the 
spectrum would still be sensitive to the CP phase and can 
be used to detect CP violation for the same hierarchy. 



In the context of a superbeam, it was recently 
pointed out [1] that the baseline of 2540 km satisfies the 
above condition for IH at a neutrino energy of 3.3 GeV 
and a narrow band neutrino beam was therefore deemed 
desirable. In this Letter, we point out that this baseline 
also satisfies the desired condition for NH, at the energy 
1.9 GeV. The two energies at which the desired condi- 
tion is satisfied are termed as magic energies, and the 
baseline is referred to as "bimagic" . The bimagic prop- 
erty, first realized in this work, makes it more desirable 
to have a broadband beam covering the range 1-4 GeV. 
We use the e-/i channel in a low energy neutrino factory 
(LENF) with a muon energy of 5 GeV 7J, as opposed to 
the ^-e channel used for superbeams 1 . The detection 
of muons is easier compared to that of electrons. More- 
over with muon charge identification, NFs do not have 
beam contamination problems, thus enabling sensitivity 
to smaller ^13 values. Thus, the bimagic nature in con- 
junction with a LENF helps in an efhcient identification 
of hierarchy, nonzero ^13 and CP violation, even with a 
single polarity of decaying muons, as we shall motivate 
and demonstrate in this Letter. It is remarkable that the 
distance 2540 km also happens to be close to the distance 
between Brookhaven and Homestake [5, as well as that 
between CERN and Pyhasalmi mine [9 , which is one of 
the proposed sites for the LENA detector. 

The bimagic baseline. — The source beam from a neu- 
trino factory that accelerates /i"*" consists of and Vf,. 
Charged current interactions at the detector can give 
muons in two ways: the original that survive as give 
/i"*" (right-sign muons) while the original v,, that oscillate 
to give /i" (wrong-sign muons). The oscillation prob- 
ability Pu^^u^ , relevant for the wrong-sign muon signal, 
can be written in the constant matter density approxi- 
mation as 1 10 
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keeping terms up to second order in a = A: 
and S13. Here Sij = sinOij, Cij = cos 9ij. Also 
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where Gi? is the Fermi constant and rig is the electron 
number density. For neutrinos, the signs of A and A are 
positive for normal hierarchy and negative for inverted 
hierarchy. A picks up an extra negative sign for anti- 
neutrinos. The last term in Eq. ([I]) clearly mixes the de- 
pendence on hierarchy and i5cp, leading to a degeneracy 
between them [1], which can be overcome if one man- 
ages to have either sin(^A) = or sin[(l — A)A] — 0. 
The first condition is achieved at the magic baseline 
(L ~ 7500 km) for all and for both the hierarchies. 
The second condition, on the other hand, is sensitive 
to hierarchy. This sensitivity can be maximized if one 
has sin[(l — A) A] ~ for one of the hierarchies and 
sin[(l — ^)A] — ±1 for the other. In such a situation, 
only the 0{a^) term in Eq. ([l]) survives for the hierarchy 
for which sin[(l — A)A] = 0, making Pe/x independent 
of both iJcp and 6*13. At the same time, for the other 
hierarchy the first term in Eq. ([I]) enhances the number 
of events as well as ^13 sensitivity, and the third term 
enhances the sensitivity to Sep- 

If we demand "IH-noCP" (no sensitivity to CP phase 
in IH), these conditions imply 
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where n, m are integers, n > 0. These two conditions 
are exactly satisfied at a particular baseline and energy, 
given by 
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Note that the relevant L is independent of any oscillation 
parameters. A viable solution for these set of equations 
(with n ~ 1 and to = 1) is L w 2540 km, p = 3.2 g/cc 
and El, = Ei^ w 3.3 GeV, as was first pointed out in 
[T]. On the other hand, one may demand "NH-noCP" 
(no sensitivity to CP phase in NH), which leads to the 
conditions 
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with n, TO integers, n 7^ and to > 0. These lead to the 
same condition on L as in Eq. (4a) except for an overall 



negative sign, while E^, continues to be given by Eq. (4b I. 
These conditions are also satisfied a.t L — 2540 km (for 
n = 1 and m — 2) at E^, = Snh ~ 1-9 GeV. The magic 
energies Em and E'nh would be suitable for a neutrino 
factory with a parent muon energy of ~ 5 GeV. 

Eqs. (4a 4b) indicate that many combinations of 
n and to are possible for a given baseline. Indeed, the 
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FIG. 1: Conversion probability Pe^ for L = 2540 km. The 
bands correspond to Sep G (0, 27r). Other parameters are 
taken as Amii = 7.65 x lO"'^ eV^ |Amii| = 0.0024 eV^ 
sin^ 612 ~ 0.3 and sin^ di2 = 0.5. The red (solid) line corre- 
sponds to ^13 = 0. 



2540 km baseline also satisfies IH-noCP at -E1H2 ~ 1-3 
GeV (n = 2, TO = 2) and NH-noCP at £;nh2 ~ o".9 GeV 
(n = 2, TO = 3). However the flux at these energies would 
be small, so we do not consider these in this Letter. 

Fig. [1] shows the probability P^^ for sin^ ^13 = 0,0.01. 
In this and all other plots, we have solved the exact neu- 
trino propagation equation numerically using the Prelim- 
inary Reference Earth Model |TT]. Clearly the IH-noCP 
and NH-noCP conditions are satisfied at the energies Em 
and -Enh, respectively. At Em, the probabilities Pe^ for 
NH and IH are distinct, hence a measurement of the neu- 
trino spectrum around this energy would be a clean way 
of distinguishing between the hierarchies. The oscillatory 
nature of Pg^ for non-zero ^13 vis-a-vis the monotonic be- 
havior for = helps in the discovery of a nonzero ^13. 
Finally, the significant widths of the bands (near Em for 
NH, and near i?NH for IH) imply sensitivity to (5cp- 

The simplified forms of probabilities at the magic ener- 
gies offer insights into the CP sensitivity at this baseline. 
At Em, we have 

Pe^(IH) « ISa^sl^cl^cl:,, 
Pe^(NH) « ISa^s^^c^^c^g + gs^gs^g 

-18V2asi2Ci2S23C23Si3 cos(Jcp + 7r/4) , (6) 
while at -Bnh, we have 

Pe^(NH) « 50ahi2cl2cl3 , 
P,^(IH) « 50ahl2cl^cl^ + {25/9)sl3sl^ 

-(50^/3)asi2Ci2S23C23Si3COs((5cp +7r/4) . (7) 

Near the magic energies, where the CP sensitivity is the 
highest, the i5cp values giving the highest and the lowest 
probabilities would be 37r/4 and 77r/4, respectively. 

Experimental setup and numerical simulation. — We 
use a magnetized totally active scintillator detector 
(TASD) which is generally used in the context of a LENF 
7]. We use a 25 kt detector with a energy threshold of 1 
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FIG. 2: Top panel: energy spectra of wrong-sign muons for 
NH (left) and IH(right). Bottom panel: number of events for 
1 year run, in the bins AEm and Ai?NH as a function of O13. 
The bands correspond to 5% error in Amli. The bands in the 
top panels also include a 10% error around sin^ 613 = 0.01. 

GeV. We choose a typical Neutrino factory setup with 5 
GeV parent muon energy and 5 x 10^^ useful niuon de- 
cays per year |12l 114] . We consider the running with only 
one polarity yCt"*" of the parent muon, so that we have a 
neutrino flux consisting of and Vg,- We assume a muon 
detection efficiency of 94% for energies above 1 GeV, 10% 
energy resolution for the whole energy range up to 5 GeV 
and a background level of 10~^ for the i^e and 
— > channels. Detection of or is not consid- 
ered in this study, which seems to have a very small effect 
when the initial flux is as large as above [H]. A 2.5% 
normalization error and 0.01% calibration error, both for 
signal and background, have also been taken into account 
throughout this study. The detector characteristics have 
been simulated by GLoBES [T5] . 

The top panel in Fig. [5] shows the energy spectra of 
wrong-sign muon events. For illustration, in addition to 
(5cp = 0, TT, we choose i5cp — 37r/4, 77r/4 which would give 
the maximum 5cv dependence near the magic energies, 
as indicated by Eqs. (|6]) and ([T]). It is clear from this 
figure that there is considerable sensitivity to CP phase 
near Em^Enu) ~ 3.3(1.9) GeV for NH (IH). It may be 
noted from the figure that the CP sensitivity for IH is 
actually better at slightly higher energies than -Enh- This 
is because the v,, spectrum at the source as well as the 
cross section of at the detector are strongly increasing 
functions of energy around E^, ~ 2 GeV, and push the 
peak in the IH spectrum to higher energies. 

In order to illustrate the effectiveness of the magic en- 
ergies, we show in the bottom panel of Fig. [2] the total 
number of events in two bins near the magic energies 
- A£;iH (3.0-3.6 GeV) and I^E^^ (1.7-2.1 GeV) - as 
functions of 6*13. Clearly, the bin Ai^in itself is enough 
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FIG. 3: The 3ct hierarchy sensitivity contours, obtained with 
a flux of 5 X 10'^^ positive muons/year on a 25 kt TASD for 
2.5 years, for true hierarchies as indicated. For parameters to 
the right of the contours, hierarchy can be identified. 



to identify the hierarchy as long as sin^ 9iz > 10^^. If the 
actual hierarchy is NH, this bin is also sensitive to 6cp- 
The sensitivity to 6*13 may be estimated by comparing 
the error bars at different di^ values. If the actual hier- 
archy is IH, one needs the events data from the energy 
bin AEjqn in order to discern Sqp and 6113. The actual 
identification of hierarchy and the measurement of Sqp 
and ^13 is done by using the complete wrong-sign events 
spectrum as well as the right-sign events spectrum. We 
present below the results of this analysis. 

Mass hierarchy determination. — In Fig.[3]we quantify 
the hierarchy sensitivity of the bimagic neutrino factory 
setup. The experimental data are generated with the 
chosen true hierarchy. The true values of sin^ ^13 and 
Sep are plotted along the axes while the true values of the 
other parameters are set to the values quoted in Fig. [T] 
For each pair of sin^ 0i3(true)-5cp(true), we obtain Xmin 
by marginalizing over other parameters. We have taken 
4% error on each of Ato2i and 9i2, and 5% error on each 
of 623 and Arri^i, for calculating the priors. 6cp has been 
varied over (0, 27r). A 2% error has also been considered 
on the earth matter profile and marginalized over. 

The contours in Fig.|3]suggest that if the true hierarchy 
is NH, then for favorable values of Sep, an exposure of 
« 3 X lO^'^ muons xkt may determine the hierarchy at 
3(7 even for sin^ ^13 as small as ^ 3 x 10"^. If the true 
hierarchy is IH then that can be established at 3a for 
sin^ ^13 > 3 X 10^*. This sensitivity is better than that 
indicated by the superbeam studies at this baseline [3 [5] . 

An optimized LENF setup with a baseline of 1300 km 
has been recently proposed [M]. However, the relatively 
small baseline does not allow matter effects to develop 
sufficiently, and one does not have the advantage of the 
magic energies. So the sensitivity of this setup to the 
hierarchy is rather limited. Indeed if the true hierarchy 
is NH, the bimagic baseline will rule out IH at 3cr for 
sin^ 9i3 values almost an order of magnitude smaller than 
the expected reach of the 1300 km setup with the same 
exposure. If IH is the true hierarchy, the performance of 
both the setups is almost the same. Thus the bimagic 
baseline is a more optimal setup as far as the hierarchy 
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FIG. 4: The 3a discovery contours for 6'i3 (upper panel) and 
CP violating phase Sep (lower panel), obtained with a flux of 
5 X lO^'^ positive muons/year on a 25 kt TASD for 2.5 years, 
for true hierarchies as indicated. For parameters to the right 
of the contours, the discovery of the relevant parameter is 
possible. 



is concerned. 

013 and Scp measurement. — The top panel of Fig. [4] 
shows that the exposure of sa 3 x 10^"^ muonsxkt will be 
able to discover a nonzero 6*13 to Scr as long as sin^ 613 > 
10~^ for either hierarchy and for any 6cp value. For NH 
and 5cp ~ 37r/4, the discovery of ^13 is possible even for 
sin^ 613 as low as 3 x 10^'"'. 

The bottom panel of Fig. |4] shows the Sep discovery 
reach with this setup. It shows that the exposure allows 
the discovery of nonzero ^cp for NH for sin^ ^13 as low 
as 10"'*, as long as 6cp ~ 37r/4. This is the Sep value 
at which we expect the highest deviation in the events 



spectrum from i5cp = 0, as indicated by Eqs. ^ and ([7|. 
For IH, the results are about one order of magnitude 
worse than those for NH. 

The discovery potential for ^13 and Sep at the bimagic 
baseline is comparable to that of the 1300 km setup if the 
true hierarchy is NH, while it is not as good if the true 
hierarchy is IH. However, note that this is valid if only /x"*" 
are available at the source. With both polarities avail- 
able, the bimagic baseline would be almost as good as 
the 1300 km setup for ^13 and Sqp, and will have a better 
sensitivity to the hierarchy. Indeed, once the hierarchy 
is identified - for which the bimagic baseline performs 
better - running the bimagic setup with /i"*" (/i~) as the 
source beam for NH (IH) would offer a sensitivity similar 
to the 1300 km setup. Thus, overall the bimagic baseline 
seems like an optimal one to probe the three most impor- 
tant unknown parameters of the leptonic mixing matrix: 
^13, Sqp and the sign of Amg^. 

Conclusion. — We have shown the "bimagic" nature 
of the 2540 km baseline: at this baseline with judi- 
cious choice of energies, the dependence of the wrong- 
sign muon signal on 613 and Sqp can be made to vanish 
for either hierarchy. This energy turns out to be around 
3.3 GeV for IH and 1.9 GeV for NH. This helps in an 
efficient identification of hierarchy even at very low 6*13, 
when one uses a neutrino factory with parent muon en- 
ergy ~ 5 GeV as a source. On the other hand the sensi- 
tivity to ^13 and ^cp is maximum at ~ 3.3 GeV for NH 
and 1.9 GeV for IH, allowing the determination of these 
parameters as well with the same beam-baseline setup. 
To exploit these features, a broadband beam of a neu- 
trino factory is more effective as compared to a narrow 
band beam. 
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